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Abstract. We investigate the sample of 1 175 new nonmagnetic DA white dwarfs with the effective temperatures > 12000 
K, which were extracted from the Data Release 1 of the Sloan Digital Sky Survey. We determined masses, radii, and bolometric 
luminosities of stars in the sample. The above parameters were derived from the effective temperatures T eff and surface gravities 
logg published in the DR1, and the new theoretical M — R relations for carbon - core and oxygen - core white dwarfs. Mass 
distribution of white dwarfs in this sample exhibits the peak at M = 0.562M Q (carbon core stars), and the tail towards higher 
masses. Both the shape of the mass distribution function and the empirical mass - radius relation are practically identical for 
white dwarfs with either pure carbon or pure oxygen cores. 
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1. Introduction 

The Sloan Digital Sky Survey is the very extensive and am- 
bitious ground-based research project, aimed at the determi- 
nation of positions, brightnesses, and acquiring of the optical 
spectra of millions of various celestial bodies. The SDSS was 
described in detail on their web page, http://www.sdss.org/, and 
in York et al. (2000). While this project was essentially de- 
signed to study galaxies, quasars, and large-scale structures in 
the Universe, its secondary targets are stars of various spectral 
types. We take advantage of the SDSS data on white dwarf stars 
of various classes, which recently were made publicly available 
as the subset of the Data Release 1 . 

Comprehensive analysis of the photometric and spec- 
tral data on white dwarfs, which are available in the 
DR1, was presented in Kleinman et al. (2004). They 
have found 2935 white dwarfs of various types and re- 
lated stars, and only few of them were spectroscopically 
identified previously. The number of well identified white 
dwarfs in their sample approaches the number of white 
dwarf stars known previously, see McCook & Sion (1999), 
and the newest version of their catalogue available at 
http://www.astronomy.villanova.edu/WDCatalog/index.html/ 

Kleinman et al. (2004) performed fitting of the observed 
spectra of 1833 nonmagnetic DA white dwarfs to the widely 
used theoretical grid of pure hydrogen model atmospheres by 
Detlev Koester (see Finley, Koester & Basri 1997). They have 
obtained effective temperatures in the range 7220 < T e ff < 
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93855 K. Surface gravities of stars in the sample range from 
low values corresponding to subdwarfs to the most massive 
white dwarfs of the highest logg = 9.0. 

The referee pointed out, that the SDSS DR1 values of logg 
determined for stars with T e ff < 12000 K are systematically 
larger than for hotter stars, and are larger than expected for 
white dwarfs, cf. Fig. 9 in Kleinman et al. (2004). The explana- 
tion of this effect was not found, however, they favor the possi- 
bility that cooler model atmospheres exhibit systematic errors 
in the parametrisation of convection, or in the fitting procedure. 
As the result, they do not claim that the high values of log g are 
real. To avoid distortion of the high-mass tail presented in Fig. 1 
we analyse a subset of the original SDSS DR1 sample which 
includes 1 175 DA white dwarfs with r eff > 12000 K. 

The Koester's grid does not include spectra of models with 
high gravities, logg > 9.0. Therefore, available determinations 
of logg also do not exceed 9.0, and this is an artificial bias 
which influences the determination of log g for over 20 most 
massive DA white dwarfs in the whole SDSS DR1 sample. 

One can note, that white dwarfs of the highest gravity 
exceeding 9.0 appear among previously known nonmagnetic 
DA stars. Compilation of the existing mass determinations by 
Nalezyty & Madej (2004) lists 6 such white dwarfs, with the 
maximum logg = 9.12. Also the most recent paper by Dahn 
et al. (2004) reports the measurement of logg = 9.46 in the 
apparently most massive star LHS 4033. Inferred mass of this 
nonmagnetic DA star is M « 1.33 M , and the corresponding 
small radius R - 0.00368 R & , place this star close to the limit- 
ing Chandrasekhar mass. 
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Fig. 1. Mass distribution of the 1175 DA white dwarfs of the 
SDSS DR1 sample with T eS > 12000 K. The peak mass is 
located at M = O.562M . Individual stars were grouped in bins 
of 0.025M o width. 
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Fig. 2. Mass - radius relation for 1175 hot white dwarfs of 
the same sample. The relation does not extend to masses M 
higher than 1 .2 M , or rather to the highest surface gravities 
logs > 9.0. 



2. Mass distribution of DA white dwarfs 

Masses M and radii R of the 1833 DA white dwarfs in the 
whole SDSS DR1 sample were directly determined from the 
corresponding values of r er y and logg. This method was in- 
troduced by Bergeron, Saffer & Liebert (1992), and is recom- 
mended for the analysis of large samples of stars. Resulting 
spectroscopic masses and radii, and their errors, strongly de- 
pend on the accuracy and reliability of the stellar atmosphere 
and the hydrogen line formation theory. 

Given r eff and log g of a star, the mass M and radius R were 
computed from the combination of the relations: 

1. Definition of the surface gravity g 

R 2 

M= — 10 log «. (1) 
G 

2. Theoretical mass - radius relation 

M = M(R, T eB ) . (2) 

We have used here M - R relations by Panei et al. (2000), 
which are briefly discussed in the following Section. The ref- 
erence set of spectroscopic M and R values was obtained as- 
suming M - R relations for pure carbon core and hydrogen 
envelope of all white dwarfs. Assumed mass of the envelope, 
M H /M = 10~ 5 . 

Computational results of our research are displayed in Figs. 
1-2. Fig. 1 presents the mass distribution of 1175 DA white 
dwarfs with T e g > 12000 K. Stars are grouped in relatively 
narrow bins of 0.025 M Q width. Masses of investigated stars 
range from M = 0.2 M Q up to 1 .2 M e . The mass distribution ex- 
hibit both the peak at 0.562 M and the high mass tail. Distinct 
asymmetry of the mass distribution and the apparent excess of 
stars just below the maximum mass of 0.562 M Q seems a real 
feature, due to the very high number of white dwarfs in Fig. 1. 



The secondary peak at 1 .2 M obviously is not a real fea- 
ture. Due to limitations in Koester's grid of model atmospheres, 
Kleinman et al. (2004) have arbitrarily assigned the maxi- 
mum logg = 9.0 to all white dwarfs, which probably have 
even higher surface gravities. Stars with logg = 9.0 clump at 
ss 1 .2 M© in the M -R relations for carbon - core white dwarfs 
(Panei et al. 2000). 

Fig. 2 presents the empirical mass-radius relation for the 
same set of 1175 hot DA stars in the SDSS DR1 sample. The 
relation does not extend to masses higher than 1 .2 M Q , or rather 
gravities logg > 9.0. For the lowest masses we note a scatter of 
the points. This is because here the surface gravities logg are 
much lower than in the high mass tail of white dwarfs, hence 
the thickness of atmospheres (contributing to R) is much larger, 
and is strongly dependent on the T e g of a star. 

3. Theoretical M-R relations 

The mass - radius relations used in this work are those de- 
scribed in Panei et al. (2000). These authors have presented 
detailed finite temperature mass - radius relations for white 
dwarfs with helium, carbon, oxygen, silicon, and iron cores 
based on a complete and updated treatment of the white dwarf 
evolution. Constitutive physics of the stellar models includes 
the equation of state for the low-density regime from Saumon 
et al. (1995) for hydrogen and helium plasmas, and ionic, 
coulomb interactions, partially degenerate electrons, electron 
exchange, and Thomas-Fermi contributions at finite tempera- 
ture for the high density regime. Radiative opacities are from 
OPAL (Iglesias & Rogers 1993) for zero metallicity and con- 
ductive opacities and neutrino emission rates for helium, car- 
bon, oxygen, silicon, and iron composition from Itoh and col- 
laborators (see Althaus & Benvenuto 1997). 

The set of evolutionary sequences presented in Panei et 
al. (2000) are particularly appropriate for our purpose because 
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they constitute a homogeneous sequence of mass-radius rela- 
tions. In particular, such relations have been constructed for 
white dwarf masses ranging from 0.15 M Q to 0.5 M Q at in- 
tervals of 5 % for helium core composition; from 0.45 M Q to 
1.2 M at intervals of 0.01 M for carbon, oxygen and silicon 
cores; and from 0.45 M to 1.0 M at intervals of 0.01 M G for 
the case of an iron core. 

In particular, for very low mass white dwarfs we have as- 
sumed a helium core composition because their progenitors 
never attain high enough central temperatures for helium to be 
ignited at their cores. The evolutionary sequences were com- 
puted down to log L/L Q = -5, and cover the effective tempera- 
ture range from 145000 K to 4000 K. In addition, two cases for 
the hydrogen envelope have been considered: M H /M = 10" 5 
(M H /M = 3 x 10" 4 for the helium core models), and M H /M = 
0. The mass of the pure helium layer above the core amounts 
to 1 % of the stellar mass. 



4. Mass and radius errors 

The SDSS DR1 tables list errors of their T e g and \ogg values 
for each star, which we denote here as ST e ff and Slogg. We 
transform them to mass and radius errors, SM and SR, in the 
following way. 

Eqs. 1 -2 imply, that small increments of all 4 parameters 
can be expressed as 



SM = 



2R 10 log « 



■SR + 



fl 2 10 log «lnl0 



SM ■■ 



dM dM 

— SR + ST eS . 

dR dT eB 



(3) 



(4) 



Eqs. |3J-|4| form the system of two algebraic equations for 
two unknown numbers, SM and 5R. Errors SI \s and 5\ogg 
are known from the input catalogue by Kleinman et al. (2004), 
and all the coefficients (including both partial derivatives) can 
easily be computed at the points (M,R), and the latter were 
determined already. 

Solving of Eqs.0-|3]yields the explicit expressions for both 
errors SR and SM (cgs units) 



(d2 a \/t 
— 10 logs lnl0£logg- ST eS 
G BS dT eff 



dR G 



d2 

SM = — 10 log «lnl0 5\ogg+ — 10 log «cW. 
G G 



(5) 



(6) 



Finally, we accepted absolute values of the above expressions 
as the real SM and SR errors. 

The average mass and radius errors resulting from our anal- 
ysis of the SDSS DR1 catalogue of DA white dwarfs are: 
(SM) = O.O59M , and (SR) = 792 km, respectively. 
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Fig. 3. The difference between DA white dwarf masses ob- 
tained from M - R relations corresponding either to pure car- 
bon, or pure oxygen cores. All the 1175 stars were grouped 
in bins of the width 0.025M o . Dots correspond to the average 
difference M(C core) - M(0 core) in each group. 



5. Core composition of white dwarf stars 

White dwarfs are thought to have a carbon/oxygen core, the 
exact proportions of which are uncertain and fixed by the mix- 
ing processes occurred prior to the formation of a white dwarf. 
Unfortunately, theoretical M-R relations for white dwarfs hav- 
ing cores with mixtures of carbon and oxygen were not avail- 
able at the time when our research was performed. However, 
we believe that the pure carbon and pure oxygen core composi- 
tions of Panei et al. (2000) would bracket the real one. Also, for 
stellar masses larger than 1 .05 solar masses, the core composi- 
tion is probably a mixture of oxygen and neon, being oxygen 
the dominant one, and for M lower than 0.45 M a pure helium 
core is expected. 

We have determined a similar set of white dwarf masses 
and radii assuming M-R relations corresponding to the pure 
oxygen core plus H envelope of the same mass as for the pure 
carbon core, Mn/M = 10~ 5 . Resulting differences between 
both grids of masses are very small. The mass of a carbon core 
DA white dwarf was always higher than the oxygen core mass 
of the same star in the whole range of 0.2 -1.2 M . The average 
difference between masses of both types is equal to 0.00592 M 
in our sample of 1 175 stars, see Fig. 3. 

All the stars in our sample were grouped in bins of the 
width 0.025 M , similarly as in Fig. 1. Dots in Fig. 3 corre- 
spond to the average difference M(C core) - M(0 core) in each 
bin. One can note, that there is little evidence of a systematic 
trend in the run of the difference M(C core) - M(0 core) vs. 
mass M(C core). The above difference increases with a white 
dwarf mass only insignificantly. 

Similar comparison was also performed for DA white 
dwarfs with carbon masses M < 0.45 M . We have found 71 
such stars in our sample of 1 175, and computed the differences 
M(He core) - M(C core) between masses derived either from 
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the M - R corresponding either to pure helium, or pure carbon 
cores , respectively. Spectroscopic mass of a helium core white 
dwarf was always higher than the mass of a carbon core star, 
and the average difference amounts to 0.0376 M , 

6. Concluding remarks 

1. We have determined the distribution of DA white dwarf 
masses for the magnitude limited sample of 1175 stars in the 
SDSS DR1 table. Our mass distribution exhibits the peak at 
M = 0.562 M Q , and slight asymmetry towards both higher and 
lower masses. 

The peak mass of this paper is in excellent agreement 
with earlier studies, which used much less numerous sam- 
ples and various methods of white dwarf mass determination. 
We present here a brief listing of previously determined peak 
masses of DA white dwarfs. Unfortunately, the early paper by 
Shipman (1979) does not provide data corresponding to data of 
other papers, which are collected in Table 1 . 

There exists a striking consistency of all peak mass de- 
terminations, in spite of relatively nonnumerous samples in 
the above listing. Our research adds the peak mass M pea k = 
0.562 M G to the list (carbon core stars), which perfectly con- 
firms the value of M pea k by Bergeron, Saffer & Liebert (1992). 
Moreover, we claim that have reliably determined the shape of 
mass distribution of DA white dwarfs, except for the single bin 
in Fig. 1, corresponding to the most massive stars. 

2. We note, that spectroscopic masses M and radii R in the 
present and earlier papers can still exhibit systematic errors 
due to uncertainties in the present theories of plasma physics. 
Theoretical profiles of the Balmer lines in spectra of DA white 
dwarfs are very strongly broadened, and the dominant contrib- 
utor is the pressure broadening of line wings. The correct de- 
scription of the pressure broadening in dense plasma of high 
gravity white dwarf atmospheres is a very difficult problem. 

Recently Wujec et al. (2002) proposed the new method 
of Stark broadening calculations based on simulation tech- 
niques. Their theory is able to compute the broadening, shift, 
and asymmetry parameters of Stark broadened hydrogen line 
opacity profile. Previous Stark broadening calculations were 
not able to include all the three plasma effects simultaneously. 
However, no grids of theoretical Balmer line profiles in white 
dwarf spectra by Wujec et al. (2002) are currently available. 

3. Catalogues of masses, radii, resulting bolometric luminosi- 
ties for both C - core and O - core white dwarfs, and all 
the input parameters are available from authors on request. 
Bolometric luminosities and their errors are given by the fol- 
lowing expressions 



Lboi = AnR 2 cr md T 



eft ' 



6L boi = SnR cr rad r e 3 ff (r eff 6R + 2R 6T eS ) 



(7) 
(8) 
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Table 1. Peak masses of the published DA white dwarf mass 
distributions, in Mp ea k/M . 

- 0.58 Koester, Schulz & Weidemann (1979), average 

mass of 122 DA white dwarfs 

- 0.603 Weidemann & Koester (1984), 70 DA WD stars 

- 0.571 McMahan (1989), 53 DA white dwarfs 

- 0.60 Weidemann (1990) 

- 0.562 Bergeron, Saffer & Liebert (1992), 129 DA white 

dwarfs 

- 0.56 Liebert & Bergeron (1995), 200 white dwarfs from 

the Palomar Green survey (Green et al. 1986) 

- 0.570 Finley, Koester & Basri (1997), 174 DA white 

dwarfs, some with cool companions 

- 0.55 Marsh et al. (1997a,b), 89 stars of the ROSAT 

All-Sky X-ray and EUV Surveys 

- 0.56 Vennes et al. (1997), 1 10 EUV selected DA WDs 

- 0.57 Vennes (1999), 141 EUV/X-ray selected WDs 
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